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Abstract 

Using Gaussian wave packet solutions, we examine how the kinetic energy is distributed in time- 
dependent solutions of the Schrodinger equation corresponding to the cases of a free particle, a 
particle undergoing uniform acceleration, a particle in a harmonic oscillator potential, and a system 
corresponding to an unstable equilibrium. We find, for specific choices of initial parameters, that 
as much as 90% of the kinetic energy can be localized (at least conceptually) in the 'front half of 
such Gaussian wave packets, and we visualize these effects. 



1 



I. INTRODUCTION 



The study of time-dependent solutions of the Schrodinger equation for initial wave packets 
of Gaussian form has a long history, dating back to the earliest days of the development of 
quantum theory. Schrodinger himself [l| and others used such solutions to discuss the 

connections between the quantum and classical formulations of mechanics and found explicit 
wave packet solutions to many of the standard problems of classical mechanics, including 
the free-particle, case of uniform acceleration, the harmonic oscillator, and the particle in 
a uniform magnetic field. These examples first appeared in a number of textbooks 
less than a decade later. 

In contrast, many treatments of introductory quantum mechanics today focus almost 
exclusively on time-independent, energy eigenstate (or stationary state) solutions of the 
Schrodinger equation, although a large number of texts do include the explicit construction 
of a time- dependent Gaussian wave packet solution for the free-particle case. This easily 
obtained, analytic result (seemingly first obtained by Darwin 1^) illustrates, in a closed form 
solution, many of the notions central to wave packet spreading and the essentially classical 
time-development of the expectation value as {x)t = (p)ot/m + {x)o. In many modern 
software packages, these solutions are also visualized, helping to illustrate the correlated 
time-development of the real and imaginaryparts, and especially their phase relationships, 
often using a popular color-coding scheme jj, ^j, jlol |. 

An example of a static image representing such time-development is shown in Fig. 1. 
One can easily note the non-symmetric pattern of the 'wiggliness' apparent for t > in 
both the real (dotted) and imaginary (dashed) parts of ^/'(a;, t), which, as we will argue 
below, can be used as a qualitative measure of the local kinetic energy. The faster (higher 
momentum, and hence 'wigglier') components are more obviously associated with the 'front 
end' of the spreading wave packet, while the 'back end' exhibits much less rapid spatial 
variation, consistent with the slower (low momentum) components trailing behind. 

In this note, we attempt to extend and amplify upon this type of intuitive observation 
and try to answer such questions as Where is the kinetic energy localized in a wave packet? 
To this end, we begin with the familiar observation that the expectation value of the kinetic 
energy operator in any time-dependent position-space state, ■i/j{x,t), can be written, using a 
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simple integration- by-parts (or IBP) argument, in the form 



2m 2m ' dx"^ 

IBP 



or 



(2) 



(T) = — / dx 
^ ' 2m ' 



2m dx dx 

dijj{x, t) 



oo 



dx 

where the 'boundary terms' (involving ip and dip/dx at ±oo) are assumed to vanish for any 
appropriately localized solution. 

This global identification suggests that we define a local kinetic energy density, T{x,t), 

via 



di:(x,t) ^ 



/ + 00 
T{x,t)dx = T{t) . (3) 
-oo 



dx 

This quantity, which is clearly locally real and positive-definite, can then be used to quantify 
the distribution of kinetic energy, and how it changes with time, for any time- dependent 
solution of the Schrodinger equation. For Gaussian solutions, we will be able to perform 
many of the desired integrals in closed form, leading to explicit analytic results for such 
related quantities as 

T^+\t)= / r{x,t)dx and T^-\t) = / r{x,t)dx (4) 

J {x)t J — oo 

where the expectation value {x)t serves to define the 'center' of the wave packet. These 
quantities can be respectively associated with the kinetic energy in the 'right' (or 'front', 
at least for packets moving generally to the right) half and the 'left' (or 'back') half of the 
wave packet. 

In what follows, we will present several explicit closed-form examples of the calculation of 
T'^^\t) for Gaussian wave packet solutions, starting in Sec. |n] with free-particle wave packets 
for a standard Gaussian momentum distribution, while in Sec. lIIII we illustrate similar results 
for Gaussian solutions to the problem of a particle undergoing uniform acceleration, working 
initially in momentum-space. In Sec. IV we use standard propagator techniques to examine 
Gaussian wave packet solutions for the harmonic oscillator problem in this context, while 
in Sec. V we extend these results to the case of an 'inverted oscillator', corresponding to a 
particle in unstable equilibrium. 



II. FREE-PARTICLE GAUSSIAN WAVE PACKETS 



The time-dependent Schrodinger equation for the one-dimensional free particle case can 
be written, and easily solved, in either position- or momentum-space in the equivalent forms 

-ih — -- — or —(t)[p,t) = ih — -— . (5) 



2m dx^ dt 2m dt 

Since we will find the momentum-space approach more useful in Sec. IIIII for the case of 

uniform acceleration, we will also use that approach here and write 

0(p,t) = 0o(p)e-^^'*/2"^^ (6) 

where 0) = (f)o{p) is the initial momentum distribution. Using this very general form, 
and the appropriate operator form of £ = ih{d/dp), we recall that 



{x)t = J (l)*{p,t) (^ih-^ (j){p,t)dp 



hoc ^ 



<Poip)^Mp) dp + — / p\Mp)\ dp (7) 



m 

= (4)0 + ^. 

m 

The position-space solution can be written, of course, using the Fourier transform as 

1 r+oo 

^(x, t) = / e^P^/^ (j){p, t) dp . (8) 



V2nh. 

The standard initial Gaussian momentum-space distribution, which gives arbitrary initial 
momentum (po) and position (xq) values, can be written in the form 



0) = 00 (P) = e-^'^^'-P")'/^ ^-^p.o/n (g) 

'71 



which gives 

{p)t = Po , {P^)t = Po + TTT ' and Apt = Apo = ■ (10) 

The explicit form of the position-space wave function is given by the Gaussian integral 



1 I — 7v — 



which can be evaluated in closed form (using the change of variables q = p—po and standard 
integrals) to obtain 

^(-^j f) = pipo{x-xo)/h p-ip'^t/2mh p-{x-xo-pot/mf/2{ahf{l+it/to) 

' ^ah{l + tt/to) 



where to = rnha^. (This result is sometimes attributed to Darwin P].) The corresponding 
probabihty density is easily shown to be 

1 



P{x,t) 



where 



x(t) = Xq + pot/m and /^t = a^y 1 + t^/tp 
and the time- dependent expectation values of position are 

{x)t = x{t) = xo+pot/m, {x'^)t = {x{t) f + ^ , and Axt 



(13) 



(14) 



(15) 
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all of which are familiar results. 

Turning now to the kinetic energy distribution defined in Eqn. (JHI), we find that the 
required spatial derivative is given by 

dip{x,t) f ipo {x — x(t)) 



dx 



{ahy{l + it/to] 



i){x,t) . 



(16) 



The kinetic energy density can therefore be written in the form 



2{x - x{t))po 



t/to 



+ 



{x - x{t)f 



'a^hy{i + tyti) 



{x,t)\\ (17) 



The expectation value of the kinetic energy is correctly given by 

f + CXD 



(18) 



and receives non-zero contributions from the first and last terms in brackets in Eqn. (fTTj) . 
since the middle term vanishes (when integrated over all space) for symmetry reasons. On 
the other hand, the individual values of T*^^^(t) in Eqn. ^ can also be evaluated giving 



2po 



ay 



+ 



(19) 



both of which are easily seen to be positive definite, as they must, due to the non- negativity 
of T{x, t). The time-dependent fractions of the total kinetic energy contained in the (+) / (— ) 
(right/left) halves of this archetypical wave packet are given by 



i?W(t) 



T(±)(t) 



± 



T{+)(t)+T(-)(t) 2 Vv^y \{2{poay + l)J ^1 + tV^o 



(pott) 



t/to 



(20) 
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which will clearly increase/decrease monotonically as t/to — > oo. The limiting values are 
then 

which have the extremal values 

^^i(^Ao^oo) = i + ^«.0.9 and R^J^it/t^ ^ ^) = ]- - ^ ^ Q.l (22) 



2 ■ " - ' 2 ^ 

which occur when 

Poa = or po = ^7= = Apo • (23) 

v2 CKV 2 

Thus, as much as 90% of the total kinetic energy of this Gaussian wave packet solution can 
be in the 'front half of the wave at long times (i.e. those for which t » to.) 

To illustrate this effect, we plot in the left column of Fig. 2, the modulus {\ilj{x,t)\, 
solid), and real (dotted) and imaginary (dashed) parts of a typical solution corresponding 
to po = (top), po — Apo (the extremal value, middle), and po = 4Apo (bottom) for long 
times {t — lOto). The values are plotted in terms of the variable x{t) — Xq +pot/m for 
easier comparison. For the po — case, it is clear that the larger momentum components (in 
magnitude) spread faster, but uniformly, in the opposite +x and —x directions, giving 
equal 'wiggliness' on each side, while for large values of po, the total kinetic energy is 
clearly increased (many more 'wiggles' everywhere), but the amount on each side of the 
expectation value peak at {x)t is roughly the same. For the extremal value of po = Apo, 
there is the clearest distinction between the 'front' and 'back' halves, as the magnitudes of 
the momentum components in the 'back' half are at a minimum, resulting in the greatest 
separation between the kinetic energy in the two halves. 

Wc can also compare the distribution of probability, described as usual by P{x,t) = 
\ilj{x,t)\'^ , to how the kinetic energy is localized. The kinetic energy density, T{x,t), can be 
scaled to the total (and possibly time-dependent) value of T{t) via 

S{x, t) = ' which satisfies / S{x, t) dx = 1 (24) 

^ ['') J-oo 

and so is normalized in the same way as the probability density. Wc plot both P{x^ t) (solid 
curve) and S{x, t) (dot-dashed curve) in the right hand cohimn of Fig. 2 for the same three 
cases considered above and we note how the shape of S'(a;, t) is correlated with the obvious 
'wiggliness' shown in the real and imaginary parts of ip{x,t). 
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III. UNIFORM ACCELERATION 



The problem of a particle under the influence of a constant force is a staple in classical 

mechanics, and was considered early in the history of quantum mechanics j^, but is less often 

discussed in introductory treatments of the subject, especially in terms of time-dependent 

solutions. For that reason, we briefly review the most straightforward momentum-space 

approach to this problem. In this case, where the potential is given by V{x) = —Fx, we can 

write the time-dependent Schrodinger equation in momentum-space as 

dj 
dp 



<P{p,t) = zh^^^ (25) 



or 

We note that the simple combination of derivatives guarantees that a function of the form 
$(p — Ft) will make the left-hand side vanish, so we assume a solution of the form t) = 
$(p — Ft)(j){p), with arbitrary and (f){p) to be determined. Using this form, Eqn. (j26p 
reduces to 

with the solution 

= e-'P'/^™^^ (28) 
We can then write the general solution as 

t) = $(j9-Ft)e-'^''/^"^^'^ (29) 

or, using the arbitrariness of as 

t) = 0o(p-Ft)e^«^'-^*)'-P')/6'"^'^) (30) 

where now (j)o{p) is some initial momentum distribution since (j){p, 0) = 0o(p)- Note that be- 
cause the p^ terms cancel in the exponential, we will be able to explicitly integrate Gaussian 
type initial momentum-space waveforms. 

For any general initial (f)o{p) we now have the time-dependent expectation values 

{p)t = {p)o + Ft (31) 

{p')t = {p')o + 2{p)oFt+{Ftf (32) 
^ {p)ot ^Ft^ 

{x)t = {x)o + + — (33) 

m 2m 



which also give the expectation value 



+ n^)\ = ^-F(x)o (34) 



2m I 2m 

which, in turn, also agrees with a similar calculation of {E)t using E = ih{d/dt), all of which 
are consistent with a particle undergoing uniform acceleration. 

Using the standard initial Gaussian momentum-space wavefunction in Eqn. Q as the 
(f)o{p — Ft) in Eqn. ((201), we can evaluate the position-space solution using Eqn. (jH)) to obtain 



ip{x,t) 



^iFt{xo-Ft^/Qm)/h ^i{po+Ft){x-xo-pot/2m)/h 



a/ y/nah{l + it /to) 

^ ^-{x-{xo+pot/m+Ft2/2m))2/2{ah)^{l+it/to) ^35^ 



The corresponding probability density is given by 

P(a;,t) = ^e-(^-"«)'/^? (36) 

- / N Pot Ft^ 

x{t) = xo + — + — 37 
m 2m 

{x)t = x{t), {x')t = {x{t)? + Y^ and = ^ (38) 

so that this accelerating wave packet spreads in the same manner as the free-particle Gaus- 
sian example. The calculation of the kinetic energy density proceeds exactly as in Sec. ITU 
with 

r(.,t) = ^f(p„ + R)^+2(p„ + K)(.-4W) 



where 



and 



2m V ' ' ' ha^ (l + ^V^o) 

^"^"^^^^^^ ^m^,t)\' (39) 



{a^hyii + tyti) 

The corresponding (+) and (— ) kinetic energies are then derived in the same way as before 
and are given by 

.-(.,^i(i)((..«,= .(^)(^).i) .0, 

which is the same result as in Eqn. with po —>■ Po + Ft. This similarity in form 

implies that the maximum (minimum) values of T^'^^t) are once again given by Eqn. ()22|l 
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(provided that t/tQ » 1) which now occurs when |j9o + Ft\ = ApQ. In this more dynamic 
situation, if po and F have different signs (so that the motion includes one 'back' and one 
'forth' component), this situation can occur twice during a single trajectory, with the roles 
of R^'^^t) and R^^\t) changing between the 'back' and the 'forth' traversal. 



IV. SIMPLE HARMONIC OSCILLATOR WAVE PACKETS 



The third model system which is easily shown to exhibit time-dependent Gaussian wave 
packet solutions is also the most frequently studied of all classical and quantum mechanical 
examples, namely the simple harmonic oscillator, defined by the potential energy V{x) = 



muo'^x^ /2. In this case, the initial value problem is 
Gaussian wave packets, by propagator techniques 

"+00 



Derhaps most easily solved, especially for 
11 1 - In this approach, one writes 



^{x,t) = dx'^{x',0)K{x,x']t,0) (41) 

J —00 

where the propagator can be derived in a variety of ways and can be written in the form 



K{x,x';t,0) 



muo 



27rihsm{ujt) 



exp 



imuo 



{{x^ + {x'Y) cos(cjt) - 2xx') 



(42) 



For the initial state wavefunction we will use position-space version of Eqn. ©, but for 
notational and visualization simplicity, we will specialize to the case of = 0, namely 

1 



V'(a;',0) 



(43) 



where (3 = ah and Axq = P/\/2. In this state, the expectation value of the energy is 



(44) 



The integral in Eqn. ()4H] can be done in closed form for the initial Gaussian in Eqn. ()4H|I 
with the result 

imujx'^ cos{ujt)'\ 1 r imAjj(3 {x — Xs{t)Y 



ijj{x, t) = exp 



2hsm{ujt) 



exp 



where 



A{t) = pcos{ujt) +i 



h 

muj[3 



sin (cut) and Xgit) = 



2hsm{ujt) A{t) 
Po sin(co't) 



mu 



(45) 



(46) 
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(We note that in the force-free hmit, when u ^ 0, one can show that this solution reduces 
to the free-particle form in Eqn. ()12|). with xq = 0, as expected.) The time- dependent 
probability density is given by 

P{x,t) = mx,t)\' = -^^l^e-(— (47) 

where 

\A{t)\ = ^Jf3^ cos2(cut) + {h/mujpy sm\ujt) (48) 

and 

{x)t = Xs{t) and Ax* = . (49) 

V2 

So, while the expectation value of the wave packet oscillates in a way which mimics the 



classical result, the time-dependent spatial width of the packet changes in time, consistent 
with ve^ gene.a, expectatio,. fo. the oscmato. case Q, [1 Th| ™ general ti„te- 
dependent Gaussian solution is described in several textbooks 21, UM and special cases of 
it are often rediscovered |2fll |. 

We note that in the special case of the minimum uncertainty wavepacket where 



P'=(^X or p=J^^p, (50) 
\mujp J V muj 

the time-dependent width of the Gaussian simplifies to 

Axt = ^ = Axo (51) 

so the wave packet oscillates with no change in shape. This is the special result seen more 
standardly 21[ - 27[ in textbooks and pedagogical articles, and is similar to the coherent- 
state like solution discussed by Schrodinger |l| in a famous paper. 

The time-dependent expectation value of the potential energy is easily found to be 



XV 1 9/ 9v 1/^9 \ . 9/ X mu'^ B"^ cos^ (cut) 
{Vix))t = -mu'{x')t = ^ (^Po + ^ ) sm^M) + . ■ (52) 



The time-dependent (total) kinetic energy then follows directly from this equation combined 
with Eqn. ()44|1 and is given by 

^, , 1,9, 1/9 \ 9/ X mu'^ B"^ sw? iut) , . 

m = 7^{V% = TT- + 77^ cos^M) + \ ' ' . (53) 



2m' ' 2m \ " 2(3 
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In order to determine the distribution of kinetic energy, however, we evaluate T{x, t) using 

dl/j{x,t) ( iVTLLO 



dx 



hsai{ijjt) 



a;COsM) - ^f^^ ) ^(X,t) 



Ait) 



and we find that 



T^^\t) 



T{t) ^pQUsm.[ijjt)cos'^{ujt) 



(54) 



(55) 



2 2^\A{t)\ 

The asymmetry {T^^\t) 7^ T'^~'\t)) in the kinetic energy distribution vanishes at half- 
integral multiples of the classical period r = 27r/c(j, namely when t = nr/2, but also at the 
classical turning points, i.e., at odd multiples of r/4. There is also no asymmetry in the case 
when pq = and the wave packet expands and contracts uniformly in both the +x and —x 
directions. Finally, there is no asymmetry in the special case of the 'fixed width' Gaussian, 
when (3 = (3o, and this property can perhaps help explain some of the special features of 
that minimum-uncertainty state. 

To better visualize the general result in Eqn. ()55p . we plot in Figs. 3 and 4 representations 
of both the wave function (modulus, real, and imaginary parts) as well as the probability 
(P(x, t)) and (scaled) kinetic energy {S{x,t)) distributions over the first quarter period. We 
note that for wave packets initially moving to the right (po > 0) as shown here, narrow 
packets, i.e., ones with (3 < (3q, typically have more kinetic energy in the 'front' half of 
the packet (Fig. 3, middle panels), while initially wider packets have the opposite behavior 
(Fig. 4) consistent with Eqn. ()55|). 

While the time-dependence of i?*^^-'(t) (defined in Eqn. pO|l ) is more varied than for the 
simpler, non-periodic, cases considered so far, as a specific example, we can examine the 
distribution of kinetic energy at times such that 



cos(co'T) = sin(ti;T) = — ^ 



In that case we find 



for instance, when 



(po/muj) 



This has extremal values of 



(2(po/ma;)2 + (/3o7/32 + /?2 




(/3oV/32) - (3' 



^ V 27r 



when 



2 ' 2/32 

which are obvious generalizations of Eqns. ()22|) and 



(56) 



(57) 



(58) 



(59) 
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V. 'INVERTED' OSCILLATOR WAVE PACKETS FOR UNSTABLE EQUILIB- 
RIUM 



The final case we present for wliich time-dependent Gaussian wave packet solutions are 
easily obtained is a generalization of the harmonic oscillator which has been described as the 

n 

'inverted oscillator' |14,] or the 'unstable particle' l28|, corresponding to the classical motion 
of a particle at the top of potential hill, given by 

V{x) = —mu^x'^ . (60) 

Many of the results obtained using propagator techniques can be easily carried over to this 
problem with the simple identifications 

— > icD , sin((x;)f:) — ^ isinh(cl;t) , and cos(ci;t) — > cosh(co't) . (61) 

For example, the position-space probability density corresponding to the initial state in 
Eqn. P3|l is given by 

Vn\B{t)\ 

where 



B{t)\ = Jp^cosh^ut) + (h/mupy sinh^ut) (63) 



and 



{x), = Ut) = '-^^^ and Ax, = \^ (64) 

and the probability density exhibits a 'runaway' (exponential) behavior. In this case, the 
long time behavior of the kinetic energy fractions is dictated by the limits 

cosh(tut), sinh(cDt) — ^ (65) 

which gives 

with the same extremal values as in Eqn. (j^^ . when pq satisfies the equivalent of Eqn. ()59|) 
with uj u. 
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VI. CONCLUSIONS AND DISCUSSION 



In this note we have provided a detailed analysis of Gaussian wave packet solutions 
of the Schrodinger equation in several model systems, helping to elucidate some of the 
qualitative aspects of wave packet time-development and spreading often seen in standard 
visualizations, and their relationship to the distribution of kinetic energy. We have focused 
on closed-form analytic results, all of which have been obtainable due to the special nature 
of the integrals which arise for Gaussian wave packets. An obvious and interesting extension 
would be to extend these results to free particle wave packets for arbitrary non-Gaussian 
initial momentum distributions, performing the required integrals in Eqn. (jSI) numerically 
to see just how general the results in Eqns. (j22|) and (j23|) are. For example, are there 
other initial wave packets and initial conditions for which the kinetic energy can be even 
more localized that the 90%/10% fraction seen here for Gaussians? Another straightforward 
generalization is to consider two-dimensional extensions of all four systems discussed here 
(then visualized as functions of both x and y) where closed form expressions are also possible, 
as well as the related case of the charged particle (confined to a plane) subject to a uniform 
magnetic field. 

Other studies of dynamical wave packet propagation have discussed aspects of the time- 
development of ilj{x,t) which arise from effects related to the differing behavior of various 
momentum-components. Examples have included discussions of the average speed of wave 
packet components which are transmitted or reflected from a rectangular barrier 3| (due 
to the energy dependence of probabilities of transmission and reflection) as well as the 
time-dependent shape of |0(j9, t)p for a wave packet hitting an infinite wall js^ or similar 
'bouncing' systems ,31j. The consideration of the distribution of kinetic energy distribution 
in many other such time-dependent wave packet solutions might prove useful in understand- 
ing their behavior. One such example certainly would be Gaussian-like wave packet solutions 
in the infinite square well, where additional interesting features, such as exact wave packet 
revivals, are present. 

Gaussian solutions such as those considered here can find use as examples of model 
systems in discussing other theoretical constructs, such as the Wigner quasi-probability 
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distribution 



ion 3 



42]. In that case, one defines Pw{x,p;t) via 



Pw{x,p;t) 



1 
1 



— oo 
oo 



tlj*{x + y,t)tlj{x-y,t) e^'Py/^dy 
(p*{p + q,t)<P{p-q,t) e-^^'^'l^dq 



(67) 
(68) 



which can then be evaluated exphcitly in closed form, in either position- or momentum- 
space, for all of the solutions discussed above. In an accompanying paper j3l we study the 
same four examples considered here, in yet a different context. 
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Figure Captions 



Fig. 1. Time-development of a free-particle Gaussian wave packet, illustrating the modulus 
{\ilj{x,t)\, solid), the real {Re{ip{x,t)), dotted) and imaginary {Im{ip{x,t)), dashed) 
parts, as functions of time. 

Fig. 2. The figures in the left column show the modulus (solid), and the real (dotted) and 
imaginary (dashed) parts of Gaussian wave packets at long times {t = lOto) plotted 
versus xq +pot/m, illustrating the large fraction of kinetic energy (note the pattern 
of 'wiggliness') in the front of the wave for the Pq = Apo case, as in Eqn. (j22I)- In the 
right column, we plot the probability density {\ip{x,t)\^, solid) and the (scaled) kinetic 
energy density {S{x,t) = T{x,t)/T(t), dashed) from Eqn. (j2H). 

Fig. 3. The figures in the top row show the modulus (solid), and the real (dotted) and imag- 
inary (dashed) parts of Gaussian wave packet solutions of the harmonic oscillator 
corresponding to xq = and po > with (3 = Pq/2 < Pq. Times over the first quarter 
period are shown. In the bottom row, we plot the probability density {{ipix, t)P, solid) 
and the (scaled) kinetic energy density {S{x, t) = T{x, t)/T{t), dashed) from Eqn. 
at the corresponding times. 

Fig. 4. Same as Fig. 3, but for an initial Gaussian wave packet solution of the harmonic 
oscillator with P = 2Pq > Pq. 
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